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Abstract. 

The general relativistic non-linear dynamics of a self-gravitating collisionless fluid 
with vanishing vorticity is studied in synchronous and comoving - i.e. Lagrangian 
- coordinates. Writing the equations in terms of the metric tensor of the spatial 
sections orthogonal to the fluid flow allows an unambiguous expansion in inverse 
powers of the speed of light. The Newtonian and post-Newtonian approximations 
are derived in Lagrangian form; the non-linear evolution of the system on super- 
horizon scales, leading to the so-called "silent universe", is also briefly discussed. 
A general formula for the gravitational waves generated by the non-linear evolution 
of cosmological perturbations is given: a stochastic gravitational-wave background 
is shown to be produced by non-linear cosmic structures, with present-day closure 
density ~ 10"^ - 10"^ on the scale of 1 - 10 Mpc. 



1. Introduction 

In these notes I review some of tlie material wliicli I presented in my lectures. These 
actually covered two different topics: the first was inflation in the early universe, mostly 
from a kinematical point of view, with some emphasis on the evolution of irregularities 
outside the Hubble radius, while the second was the dynamics of self-gravitating irrota- 
tional dust in Lagrangian coordinates. A strong link between the two topics is given by 
the use of General Relativity in dealing with cosmological problems on both very large 
(super-horizon) and small scales, where it is usually thought that Newtonian gravity 
should provide a detailed description of the system. Given the many excellent reviews 
on cosmological inflation, among which the lectures given by Rocky Kolb at this school 
(this volume), I preferred to devote these notes to the second topic. Most of the material 
presented here is based on a recent work in collaboration with David Terranova [0, as 
well as on many ideas which come out from a research program on "silent universes", 
in collaboration with Ornella Pantano, Diego Saez and Marco Bruni !§]• 
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The origin of my interest, however, goes back to my first reading of the lectures by 
George Elhs at the 1969 Varenna School |^, whose title, "Relativistic Cosmology", I 
borrowed for these notes. Let me finally mention that the material in these notes is 
strongly related to that in the lectures by Francois Bouchet and Thomas Buchert (this 
volume), which deal with various aspects of the dynamics of self-gravitating dust within 
the Newtonian approximation. 

The first aim of these notes is to try to convince the skeptical reader that a general 
relativistic approach is relevant for the structure formation problem in cosmology, i.e. 
that there exists a class of problems in which non-linearity and relativistic effects both 
come into play. The gravitational instability of cold matter in a cosmological framework 
is in fact usually studied within the Newtonian approximation. This basically consists 
in neglecting terms higher than the first in metric perturbations around a matter- 
dominated Friedmann-Robertson- Walker (FRW) background, while keeping non-linear 
density and velocity perturbations. This approximation is usually thought to produce 
accurate results in a wide spectrum of cosmological scales, namely on scales much larger 
than the Schwarzschild radius of collapsing bodies and much smaller than the Hubble 
horizon scale, where the peculiar gravitational potential yjg, divided by the square of 
the speed of light (? to obtain a dimensionless quantity, keeps much less than unity, 
while the peculiar matter fiow never becomes relativistic. To be more specific, the 
Newtonian approximation consists in perturbing only the time-time component of the 
FRW metric tensor by an amount 2ipg/c^, where (pg is related to the matter density 
fiuctuation 5 via the cosmological Poisson equation, V^.v?g(x, t) = 47rG'a^(t)^b(t)5(x, t), 
where is the background matter density and a{t) the appropriate FRW scale-factor; 
the Laplacian operator has been used here with its standard meaning of Euclidean 
space. The fiuid dynamics is then usually studied in Eulerian coordinates by accounting 
for mass conservation and using the cosmological version of the Euler equation for a 
self-gravitating pressureless fiuid, as long as the fiow is in the laminar regime, to close 
the system. To motivate the use of this "hybrid approximation", which deals with 
perturbations of the matter and the geometry at a different perturbative order, one can 
either formally expand the correct equations of GR in inverse powers of the speed of 
light (e.g. Ref. P]), or simply notice that the peculiar gravitational potential is strongly 
suppressed with respect to the matter perturbation by the square of the ratio of the 
perturbation scale A to the Hubble radius th = cH~^ {H being the Hubble constant): 

Such a simplified approach, however, already fails in producing an accurate 
description of the trajectories of relativistic particles, such as photons. Neglecting the 
relativistic perturbation of the space-space components of the metric, which in the so- 
called longitudinal gauge is just —2ipg/c^, would imply a mistake by a factor of two 
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in well-known effects such as the Sachs- Wolfe, Rees-Sciama and gravitational lensing, 
as it would be easy to see, by looking at the solution of the eikonal equation. In 
other words, the level of accuracy not only depends on the peculiar velocity of the 
matter producing the spacetime curvature, but also on the nature of the particles 
carrying the signal to the observer. Said this way, it may appear that the only 
relativistic correction required to the usual Eulerian Newtonian picture is that of writing 
the metric tensor in the revised, "weak field", form (e.g. Ref. However, as we 
are going to show, this is not the whole story. It is well-known in fact that the 
gravitational instability of aspherical perturbations (which is the generic case) leads 
to the formation of very anisotropic structures whenever pressure gradients can be 
neglected (e.g. Ref.|jlO|). Matter first flows in almost two-dimensional structures called 
pancakes, which then merge and fragment to eventually form one-dimensional filaments 
and point-like clumps. During the process of pancake formation the matter density, the 
shear and the tidal field formally become infinite along evanescent two-dimensional 
configurations corresponding to caustics; after this event a number of highly non-linear 
phenomena, such as vorticity generation by multi-streaming, merging, tidal disruption 
and fragmentation, occur. Most of the patology of the caustic formation process, such 
as the local divergence of the density, shear and tide, and the formation of multi- 
stream regions, are just an artifact of extrapolating the pressureless fluid approximation 
beyond the point at which pressure gradients and viscosity become important. In spite 
of these limitations, however, it is generally believed that the general anisotropy of the 
collapse configurations, either pancakes or filaments, is a generic feature of cosmological 
structures originated through gravitational instability, which would survive even in the 
presence of a collisional component. 

This simple observation already shows the inadequacy of the standard Newtonian 
paradigm. According to it, the lowest scale at which the approximation can be 
reasonably applied is set by the amplitude of the gravitational potential and is given 
by the Schwarzschild radius of the collapsing body, which is negligibly small for any 
relevant cosmological mass scale. What is completely missing in this criterion is the 
role of the shear which causes the presence of non-scalar contributions to the metric 
perturbations. A non-vanishing shear component is in fact an unavoidable feature of 
realistic cosmological perturbations and affects the dynamics in (at least) three ways, 
all related to non-local effects, i.e. to the interaction of a given fluid element with the 
environment. 

First, at the lowest perturbative order the shear is related to the tidal field generated 
by the surrounding material by a simple proportionality law. This sort of non-locality, 
however, is coded in the initial conditions of each fluid-element through a Coulomb-like 
interaction with arbitrarily distant matter. Because of its link with the initial data of 
each fluid element one can consider it as a local property. The later modification of 
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these shear and tidal fields is one of the consequences of the non-linear evolution. 

Second, it is related to a dynamical tidal induction: the modification of the 
environment forces the fiuid element to modify its shape and density. In Newtonian 
gravity, this is an action- at- a- distance effect, which starts to manifest itself in second- 
order perturbation theory as an inverse-Laplacian contribution to the velocity potential 
(e.g. Catelan et al. and references therein). 

Third, and most important here, a non-vanishing shear field leads to the generation 
of a traceless and divergenceless metric perturbation which can be understood as 
gravitational radiation emitted by non-linear perturbations. This contribution to the 
metric perturbations is statistically small on cosmologically interesting scales, but 
it becomes relevant whenever anisotropic (with the only exception of exactly one- 
dimensional) collapse takes place. In the Lagrangian picture considered here, such an 
effect already arises at the post-Newtonian (PN) level. 

Note that the two latter effects are only detected if one allows for non-scalar 
perturbations in physical quantities. Contrary to a widespread belief, in fact, the choice 
of scalar perturbations in the initial conditions is not enough to prevent tensor modes 
to arise beyond the linear regime in a GR treatment. Truly tensor perturbations are 
dynamically generated by the gravitational instability of initially scalar perturbations, 
independently of the initial presence of gravitational waves. 

Recently a number of different approaches to relativistic effects in the non-linear 
dynamics of cosmological perturbations have been proposed. Matarrese, Pantano and 
Saez proposed an algorithm based on neglecting the magnetic part of the Weyl tensor 
in the dynamics, obtaining strictly local fluid-flow evolution equations, i.e. the so- 
called "silent universe" . Using this formalism Bruni, Matarrese and Pantano studied 
the asymptotic behaviour of the system, both for collapse and expansion, showing, in 
particular, that this kind of local dynamics generically leads to spindle singularities for 
collapsing fluid elements. This formalism, however, cannot be applied to cosmological 
structure formation inside the horizon, where the non-local tidal induction cannot 
be neglected, i.e. the magnetic Weyl tensor if^ is non-zero, with the exception of 
highly speciflc initial conflgurations 0,|T^. Rather, it is probably related to the non- 
linear dynamics of an irrotational fluid outside the (local) horizon One possible 
application is in fact connected to the so-called Cosmic No-hair Theorem, i.e. to 
the conjecture that expanding patches of an initially inhomogeneous and anisotropic 
universe asymptotically tend to almost FRW solutions, thanks to the action of a 
cosmological constant-like term. The self-consistency of these "silent universe" models 
has been recently demonstrated by Lesame, Dunsby and Ellis |T^, extending an earlier 
analysis by Barnes and Rowlingson Lesame, Ellis and Dunsby |15] showed that 



any non-zero iJ'^ has non-vanishing divergence (although this appears only as a third- 
order effect in the amplitude of perturbations around the FRW background), reflecting 
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the fact that the shear and the electric tide generally have a different eigenframe. A 
local-tide approximation for the non-linear evolution of collisionless matter, which tries 
to overcome some limitations of the Zel'dovich approximation [|16|, has been recently 



proposed by Hui and Bertschinger . 

Following Matarrese and Terranova [|I| we review here a more "conservative" 
approach based on expanding the Einstein and continuity equations in inverse powers of 
the speed of light, which will then define a Newtonian limit and, at the next order, post- 
Newtonian corrections. What is new in this approach is the choice of the synchronous 
and comoving gauge, because of which the method can be called Lagrangian. Various 
approaches have been proposed in the literature, which are somehow related to the 
present one. A PN approximation has been followed by Futamase to describe 
the dynamics of a clumpy universe. Tomita also used non-comoving coordinates 
in a PN approach to cosmological perturbations. Shibata and Asada recently 
developed a PN approach to cosmological perturbations, but they also used non- 
comoving coordinates. Kasai |2T| analyzed the non-linear dynamics of dust in the 
synchronous and comoving gauge; his approximation methods are however largely 
different. Finally, in a series of papers, based on the Hamilton-Jacobi approach (e.g. 
Ref. p2| and references therein) a new approximation technique has been developed, 
which relies on an expansion in higher and higher gradients of an initial perturbation 
"seed". 



2. Relativistic dynamics of collisionless matter 

In this section we will derive the equations governing the evolution of an irrotational fluid 
of dust (i.e. p = u = 0) in a synchronous and comoving system of coordinates (actually 
the possibility of making these two gauge choices simultaneously is a peculiarity of 
irrotational dust, which holds at any time, i.e. also beyond the linear regime). The 
starting point will be the Einstein equations Rat — ^QatR = ^^Tab, with Rat the 
Ricci tensor, and the continuity equation T'^'!^ = for the matter stress-energy tensor 
T"* = gc^u"'u^, where g is the mass density and u"" the fluid four-velocity (normalized 
to u^Ua = — c^). The line-element reads 

rfs2 ^ ^c^dt^ + hapiq, t)dq''dq^ . (1) 

The fluid four-velocity in comoving coordinates is = (c, 0,0,0). A fundamental 
quantity of this analysis will be the velocity-gradient tensor, which is purely spatial, 

0"/3 = = ih^'Kp , (2) 

where a dot denotes partial differentiation with respect to the proper time t. The tensor 
0^3 represents the extrinsic curvature of the spatial hypersurfaces orthogonal to u"". 
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Thanks to the spacetime sphtting obtained in our frame, the 10 Einstein equations 
can be immediately divided into 4 constraints and 6 evolution equations. The time-time 
component of the Einstein equations is the so-called energy constraint, which reads 

- e"^e^„ + ^^^r = iqttGq , (3) 

where the volume-expansion scalar is just the trace of the velocity-gradient tensor, 
(^^i? is the trace of the three-dimensional Ricci curvature, ^^^i?"^, of the spatial 
hypersurfaces of constant time. 

The space-time components give the momentum constraint, 

6>"/?||a = 0,0 , (4) 

where greek indices after a comma denote partial derivatives, while after a double vertical 
bar they denote covariant derivatives in the three-space with metric hap- 

Finally, the space-space components represent the only truly evolution equations, 
i.e. those which contain second-order time derivatives of the metric tensor. They indeed 
govern the evolution of the extrinsic curvature tensor and read 

9"^ + 00"^ + = 47rG^r^ . (5) 

Taking the trace of the last equation and combining it with the energy constraint, one 
obtains the Raychaudhuri equation, 

+ 0"^0^^ + At^Gq = . (6) 

Mass conservation is provided by the equation 

q = -Qq. (7) 

Given that = ^h""'hja = d{\nh^^'^)/dt, where h = det hap. One can write the 
solution of this equation in the form 

^(q,^) = ^o(q)[Mq,^)/^o(q)]"'/'- (8) 

Here and in what follows quantities with a subscript are meant to be evaluated at 
some initial time to- 

Let us also introduce the so-called electric and magnetic parts of the Weyl tensor, 
which are both symmetric, flow-orthogonal and traceless. They read, respectively. 



and 



n = Is^p {e^.e", - e') + 00"^ - e^'.e'p + c' ( ^'^r^^ - ^5^^ (9) 

= lhp,{v'^'e%^s + V^'^'Q'.ws) . (10) 

where Ty^/^T/t^^/^e^/^T three-dimensional, completely anti-symmetric, Levi-Civita 

tensor relative to the spatial metric hap and 6°'^'^ is such that e^^^ = 1, etc... . 
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Notice that, while the definition of the electric tide E'^p is completely fixed, because 
of its well-known Newtonian limit, the magnetic tensor field has no straightforward 
Newtonian counterpart, and can be therefore defined up to arbitrary powers of the 
speed of light. The definition adopted here is the most straightforward one; it is such 
that no explicit powers of c appear, which means that its physical dimensions are 1/c 
those of E°'p. This choice can be motivated in analogy with electrodynamics, where the 
magnetic vector field is also scaled by 1/c with respect to the electric one. 



With the purpose of studying gravitational instability in a FRW background, it is 
convenient to factor out the homogeneous and isotropic solutions of the above equations, 
perform a conformal rescaling of the metric with conformal factor a{t), the scale-factor 
of FRW models, and adopt the conformal time r, defined by dr = dt/a{t). 

The line-element is then written in the form 

rfs^ = a2(r)[-cW + 7«/3(q,r)dg"V] , (11) 

where a^(r)7Q,/3(q, r) = hai3{q,t{T)). For later convenience let us fix the Lagrangian 
coordinates to have physical dimension of length and the conformal time variable r 
to have dimension of time. As a consequence the spatial metric ■ja/s is dimensionless, as 
is the scale-factor a(r) which must be determined by solving the Friedmann equations 
for a perfect fiuid of dust 

'a'\^ SvrG 22 

—eba - Kc , (12) 

a" /a'^2 



2 +Kc^ = . (13) 

a \a / 

Here primes denote differentiation with respect to the conformal time r and k represents 
the curvature parameter of FRW models, which, because of our choice of dimensions, 
cannot be normalized as usual. So, for an Einstein-de Sitter universe k = 0, but for a 
closed (open) model one simply has k > {k < 0). Let us also note that the curvature 
parameter is related to a Newtonian squared time-scale through kn = i^c^ (e.g. 
Refs.||23||,||24||);in other words k is an intrinsically PN quantity. 

By subtracting the isotropic Hubble-flow, one introduces a peculiar velocity-gradient 
tensor 

where = (c/a, 1, 1, 1). 

Thanks to the introduction of this tensor one can rewrite the Einstein's equations 
in a more cosmologically convenient form. The energy constraint becomes 

- d^'J^ + 4-^9 + c2(7^ - 6k) = IQ-nGa^Qb^ , (15) 
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where TZ°'p{'y) = ^^^R°'f^(h) is the conformal Ricci curvature of the three-space, 
i.e. that corresponding to the metric 7^/3; for the background FRW solution 7^^^ = 
(1 + f one has 72."g(7^'^^) = 2k5"^. We also introduced the density contrast 

6 = {g- Qb)/Qh- 

The momentum constraint reads 

= ■ (16) 

To avoid excessive proliferation of symbols, the double vertical bars are used here and 
in the following for covariant derivatives in the three-space with metric ja/s- 

Finally, after replacing the density from the energy constraint and subtracting the 
background contribution, the extrinsic curvature evolution equation becomes 
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47^"^-(7^+2/^)5"^ 



0.(17) 



The Raychaudhuri equation for the evolution of the peculiar volume-expansion scalar 
becomes 



^' + -'d + ^''J\ + A7rGa'0b5 = O. (It 
a 



The main advantage of this formalism is that there is only one dimensionless (tensor) 
variable in the equations, namely the spatial metric tensor 7^^, which is present with its 
partial time derivatives through -i?"^, and with its spatial gradients through the spatial 
Ricci curvature TZ'^p- The only remaining variable is the density contrast which can be 
written in the form 

5(q,r) = (1 + 5o(q))[7(q,r)/7o(q)]"'/' - 1 , (19) 

where 7 = det 7^/3. A relevant advantage of having a single tensorial variable, for our 
purposes, is that there can be no extra powers of c hidden in the definition of different 
quantities. 

Our intuitive notion of Eulerian coordinates, involving a universal absolute time and 
globally flat spatial coordinates is intimately Newtonian; nevertheless it is possible to 
construct a local coordinates system which reproduces this picture for a suitable set of 
observers. This issue has been already addressed in Refs. 0,0, where local Eulerian - 
FRW comoving - coordinates where introduced, related to the Lagrangian ones 
via the Jacobian matrix with elements 

J\{q, r)^^^5t + V^M. r), A = 1,2,3, (20) 

where ■D'^(q, r) is called deformation tensor. Each matrix element labelled by 
the Eulerian index A can be thought as a three-vector, namely a triad, defined on the 
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hypersurfaces of constant conformal time. As shown in Refs.p|,|^, they evolve according 
to 

J^l = ^\J\ , (21) 

which also follows from the condition of parallel transport of the triads relative to q 
along the world-line of the corresponding fluid element D{aJ^^) / Dt = 0. 

Our local Eulerian coordinates are such that the spatial metric takes the Euclidean 
form 6ab, i-e. 

7a/?(q, r) = SabJ\{ci, r) J^^(q, r) . (22) 
Correspondingly the matter density can be rewritten in the suggestive form 

Qici,T) = f5f,(r)(l + 5o(q))[J(q,r)/Jo(q)]"' , (23) 

where JT" = detj''^^. Note that, contrary to the Newtonian case, it is generally im- 
possible in GR to fix iTo = 1, as this would imply that the initial Lagrangian space is 
conformally flat, which is only possible if the initial perturbations vanish. 

We are now ready to deal with the linearization of the equations obtained above. 
This will be done mostly for pedagogical purposes, in that it will allow us to obtain 
a number of results which will turn out to be useful for the expansion. Apart 
from this, it can be interesting to re-obtain the classical results of linear theory in the 
comoving and synchronous gauge only in terms of the spatial metric coefficients. 

Let us then write the spatial metric tensor of the physical (i.e. perturbed) space- 
time in the form 

7a/3 = 7a/3 + Wa/B , (24) 

with 7q,/3 the spatial metric of the background space - in our case the maximally 
symmetric FRW one, ■ja/s = 7a/f^ ^ and Wa/3 a small perturbation. The only non- 
geometric quantity in our equations, namely the initial density contrast 60, can be 
assumed to be much smaller than unity. 

As usual, one can take advantage of the maximal symmetry of the background FRW 
spatial sections to classify metric perturbations as scalars, vectors and tensors. One then 
writes 

U!aP = XlalB + C|a/3 + ^(^a|/3 + ^I3\a) + 71"a/3 , (25) 

with 

r|a = vr"„ = 7rV = 0, (26) 

where a single vertical bar is used for covariant differentiation in the background 
three-space with metric 'ja/s- In the above decomposition x and ( represent scalar 
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modes, vector modes and vr"^ tensor modes (indices being raised by the contravariant 
background three-metric) . 

Before entering into the discussion of the equations for these perturbation modes, 
let us quote a result which will be also useful later. In the "i?"^ evolution equation and in 
the energy constraint the combination V^^ = 47^°^ — {71 + 2^)5°'^ and its trace appear. 
To first order in the metric perturbation one has 



(27) 



where V^(-) = (■)|''^ . Only the scalar mode x ^^id the tensor modes contribute to the 
three-dimensional Ricci curvature. 

As well known, in linear theory scalar, vector and tensor modes are independent. 
The equation of motion for the tensor modes is obtained by linearizing the traceless 
part of the evolution equation. One has 



^ " ^0^ -rr ' 
T^afS + 2 — TCal3 

a 



c2(V^ - 2K)7r. 



a(3 







(2^ 



which is the equation for the free propagation of gravitational waves in a FRW 
background in the Einstein-de Sitter case). The general solution of this equation is 
well-known (e.g. Ref.^]) and will not be reported here. 

At the linear level, in the irrotational case, the two vector modes represent gauge 
modes which can be set to zero, ^" = 0. 

The two scalar modes are linked together through the momentum constraint, which 
leads to the relation 



X = Xo + t^{C- Co) • 
The energy constraint gives 



(V^ + 3n) 



a 



C' + {A7rGa^Q,-Kc'){C-Co)-c\o 



while the evolution equation gives 



a 



C" + 2-C 

a 



(29) 



(30) 



(31) 



An evolution equation only for the scalar mode ( can be obtained by combining 
together the evolution equation and the energy constraint; it reads 



(V^ + 3k) 



C" + -C'-47rG'a2^,(C-Co 



On the other hand, linearizing the solution of the continuity equation, gives 



(32) 



(33) 
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which replaced in the previous equation gives 

6" + -6' - AnGa^gbS = . (34) 
a 

This is the well-known equation for linear density fluctuation, whose general solution 
can be found in Ref.|p3||. Once 6{t) is known, one can easily obtain ( and which 
completely solves the linear problem. 

Eq.(32) above has been obtained in whole generality; one could have used instead 
the well-known residual gauge ambiguity of the synchronous coordinates to simplify its 
form. In fact, ( is determined up to a space-dependent scalar, which would neither 
contribute to the spatial curvature, nor to the velocity-gradient tensor. For instance, 
one could flx (q so that (V^ + 3k)Co = — 2(5o, so that the ( evolution equation takes the 
same form as that for 6. 

In order to better understand the physical meaning of the two scalar modes x 
C, let us consider the simplest case of an Einstein-de Sitter background (k = 0), for 
which a(r) oc r^. By flxing the gauge so that V^Co = — 2(5o one obtains = Xo and 

Cir) = ^Xor' + Bor-\ (35) 

where the amplitude of the decaying mode is an arbitrary function of the spatial 
coordinates. Consistency with the Newtonian limit suggests xo = "-^Vo, with (po the 
initial peculiar gravitational potential, related to 6o through V'^ipo = AnGalgobSo. One 
can then write 

C(r) = -^^ot' + BoT-' . (36) 

This result clearly shows that, at the Newtonian level, the linearized metric is 
1ai3 = 5ap + C|a/3; whilc the pcrturbatiou mode x is already PN. Note that also the 
tensor modes are at least PN. 

These results also conflrm the above conclusion that in the general GR case the 
initial Lagrangian spatial metric cannot be flat, i.e. j7o 7^ 1, because of the initial 
"seed" PN metric perturbation xq- 



3. Newtonian approximation 

The Newtonian equations in Lagrangian form can be obtained from the full GR 
equations by an expansion in inverse powers of the speed of light; as a consequence 
of our gauge choice, however, no odd powers of c appear in the equations, which implies 
that the expansion parameter can be taken to be 

Let us then expand the spatial metric in a form analogous to that used in the linear 
perturbation analysis of Section 2: 

lap = lap + lAT + ^{^^ ' (37) 
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where the c dependence of the metric perturbation was made exphcit. The actual 
convergence of the series requires that the PN metric perturbation ^w^^^'* is much 
smaller than the background Newtonian metric 70/3. Let us first concentrate on the 
Newtonian metric; the properties of will be instead considered in Section 4. 

To lowest order in our expansion, the extrinsic curvature evolution equation and the 
energy constraint imply that = 'P'^(7) = 0, and recalling that k = Ktv/c^, 

= 7^"^(7) = : (38) 

in the Newtonian limit the spatial curvature identically vanishes (e.g. Ref. |]^). This 
important conclusion implies that 7^/3 can be transformed to 5ab globally, i.e. that one 
can write 

7a/3 = ^abJ^o^J^p , (39) 

with integrable Jacobian matrix coefficients. In other words, at each time r there exist 
global Eulerian coordinates such that 

x(q,r) = q+S(q,r), (40) 

where S(q, r) is called the displacement vector, and the deformation tensor becomes in 
this limit 

= (41) 
The Newtonian Lagrangian metric can therefore be written in the form 

One can rephrase the above result as follows: the Lagrangian spatial metric in 
the Newtonian limit is that of Euclidean three-space in time-dependent curvilinear 
coordinates g", defined at each time r in terms of the Eulerian ones by inverting 
Eq.(40) above. As a consequence, the Christoffel symbols involved in spatial covariant 
derivatives (which will be indicated by a single bar or by a nabla operator followed by 
greek indices) do not vanish, but the vanishing of the spatial curvature implies that 
these covariant derivatives always commute. 

Contrary to the evolution equation and the energy constraint, the Raychaudhuri 
equation and the momentum constraint contain no explicit powers of c, and therefore 
preserve their form in going to the Newtonian limit. These equations therefore determine 
the background Newtonian metric 7^/3, i.e. they govern the evolution of the displacement 
vector S. 

The Raychaudhuri equation becomes the master equation for the Newtonian 
evolution; it takes the form 

^' + -^ + ^^,r^ + ATiGa^g.iT^/^ - 1) = , (43) 



7..(q, r) = Sas ( Si + ^^^^^ ){S% + ^^^^^ ) . (42) 
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where 

= Irx^ , (44) 

and, for simplicity, 60 = was assumed (a restriction which is, however, not at all 
mandatory). We also used the residual gauge freedom of our coordinate system to set 
laisiTo) = 6af3, implying So = 1, i-e. to make Lagrangian and Eulerian coordinates 
coincide at the initial time. That this choice is indeed possible in the Newtonian limit 
can be understood from our previous linear analysis, where this is achieved by taking, 
e.g., Co = 0. 

The momentum constraint, 

K\, = ^.^ , (45) 

is actually related to the irrotationality assumption. 

Let us also notice a general property of our expression for the Lagrangian metric: 
at each time r it can be diagonalized by going to the local and instantaneous principal 
axes of the deformation tensor. Calling % the eigenvalues of the metric tensor, Sa those 
of the Jacobian and da those of the deformation tensor, one has 

7a(q, r) = Ji^iq, r) = (1 + J,(q, r))' . (46) 

The diagonal form of the metric tensor will be reconsidered in the frame of 
the Zel'dovich approximation. Beyond the mildly non-linear regime, where this 
approximation is consistently applied, diagonalizing the metric is in general, i.e. 
apart from specific initial configurations, of smaller practical use, because metric (and 
deformation) tensor, shear and tide generally have different eigenvectors. 

From this expression it becomes evident that, at shell-crossing, where some of the 
Jacobian eigenvalues go to zero, the related covariant metric eigenvalues just vanish. On 
the other hand, other quantities, like the matter density, the peculiar volume expansion 
scalar and some eigenvalues of the shear and tidal tensor will generally diverge at the 
location of the caustics (see Ref. |Q, for a discussion). This diverging behaviour makes 
the description of the system extremely involved after this event. Although dealing 
with this problem is far outside the aim of the present paper, let us just mention that a 
number of ways out are available. One can convolve the various dynamical variables by 
a suitable low-pass filter, either at the initial time, in order to postpone the occurrence 
of shell-crossing singularities, or at the time when they form, in order to smooth the 
singular behaviour; alternatively one can abandon the perfect fluid picture and resort to 
a discrete point-like particle set, which automatically eliminates the possible occurrence 
of caustics, at least for generic initial data. At this level, anyway, we prefer to take a 
conservative point of view and assume that the actual range of validity of our formalism 
is up to shell-crossing. 
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A more direct way to deal with the Lagrangian Newtonian equations is to write 
them in terms of the Jacobian matrix J^^. This approach is obviously related to the 
more usual ones in terms of the displacement vector S or in terms of the deformation 
tensor |25| , [gg , |28l , 



In order to rewrite the Raychaudhuri equation in terms of the Jacobian matrix, we 
notice that 

t?"/ = T^j^f, " - ^-.r^ , (47) 

where we have introduced the inverse Jacobian matrix 

(48) 

where Eulerian indices are raised and lowered by the Kronecker symbol. To make 
explicit the notation, we just stress that elements of will be characterized by a 
greek (i.e. Lagrangian) index subscript, while elements of the inverse matrix will 
be characterized by a greek index superscript. 

Replacing the latter identity into the Newtonian expression Eq.(43) yields 

J'lJ^" + -J-'J' = iiTGa^gbil - J-') . (49) 
We also notice that the parallel transport condition can be rewritten in the form 

n = J''aJ%' ■ (50) 
This equation, together with -(9"^ = |7"'^7'y/3' gives the general relation 

J^Aa ^ /3 = Aa3 /3 • (51) 

Replacing these relations in the momentum constraint we obtain in whole generality 

J'^ A^^P ||a + A\\a^^p = (<^ ^^'\fi ■ (52) 

On the other hand, in the Newtonian limit we have 

•^^a^P = ^^(i,a 1 (53) 

as it follows from the fact that = S^^^^. Using this commutation property it is 

easy to verify that 

r^-y = J\J'^fi^-^ ■ (54) 

Thanks to the latter relation and to the well-known matrix identity Trln J = IndetJ, 
it is straightforward to verify that the momentum constraint in the Newtonian limit 
becomes an identity. It is then clear that Eq.(51) is more fundamental than the 
momentum constraint: it plays the role of an irrotationality condition written in 
Lagrangian space. This is of course equivalent to the standard form [compare with 
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Eq.(59) in Ref.[^] e'^^"' J pj'j^^. This equation, together with the Raychaudhuri 
equation above, completely determines the Newtonian problem, in terms of either the 
Jacobian matrix, the deformation tensor or the displacement vector. 

As demonstrated above, it is always possible, in the frame of the Newtonian 
approximation, to define a global Eulerian picture. This will be the picture of the 
fluid evolution as given by an observer that, at the point x = q + S(q, r) and at the 
time r observes the fluid moving with physical peculiar three- velocity v = dS/dr. 
From the point of view of a Lagrangian observer, who is comoving with the fluid, 
the Eulerian observer, which is located at constant x, is moving with three-velocity 
(iq(x, r)/(ir = —v. 

The line-element characterizing the Newtonian approximation in the Eulerian frame 
is well-known (e.g. Ref . [p3l ) 



ds' = a'ir) 



1 + ) cV + 5ABdx^dx^ 



(55) 



with (yjg the pecuhar gravitational potential, determined by the mass distribution 
through the Eulerian Poisson equation, 

V^.yp,(x, r) = 47rG'a2(r)f?5(r)5(x, r) , (56) 

where the Laplacian V^, as well as the nabla operator V, have their standard Euclidean 
meaning. The perturbation in the time-time component of the metric tensor here comes 
from the different proper time of the Eulerian and Lagrangian observers. 

It is now crucial to realize that all the dynamical equations obtained so far, being 
entirely expressed in terms of three-tensors, keep their form in going to the Eulerian 
picture, only provided the convective time derivatives of tensors of any rank (scalars, 
vectors and tensors) are modified as follows: 

^ d ^ dS 

^ +v.V, v= — . 57 

Ur or dr 

This follows from the fact that, for the metric above, = T^^ = = 0, which 
also obviously implies that covariant derivatives with respect to reduce to partial 
ones. 

The irrotationality assumption now has the obvious consequence that we can define 
an Eulerian velocity potential through 

v(x,r) = V<l',(x,r) . (58) 

The Newtonian peculiar velocity-gradient tensor then becomes 

= ' ^'^^ 

because of which the momentum constraint gets trivially satisfied and the magnetic 
Weyl tensor becomes identically zero in the Newtonian limit. 
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We can now write the Raychaudhuri equation for the Eulerian pecuhar volume- 
expansion scalar {}, and use the Poisson equation to get, as a first spatial integral, the 
Euler equation 

V ' + V • Vv + -V = -V^g . (60) 

(X 

This can be further integrated to give the Bernoulli equation 

K+'^^V + \{y^,f = -iPg. (61) 

On the other hand, by taking gradients of the Euler equation we can obtain an 
Eulerian evolution equation for the tensor ^ab- More interesting is that this equation 
can be transported back to the Lagrangian frame to get 



+ + = -V^ir, ■ (62) 

where ^p^^'^ must be thought as a Lagrangian peculiar gravitational potential to be 
determined through the Lagrangian Poisson equation 

^\f^ = 47rGa'g,{T'^' - 1) . (63) 

These two Lagrangian expressions will turn out to be very useful for the PN calcu- 
lations of Section 4. 

Having shown the equivalence of this method, in the Newtonian limit, with the 
standard one, it is now trivial to recover the Zel'dovich approximation [Q. This is 
obtained by expanding Eq.(49) and Eq.(51) to first order in the displacement vector. 
The result is 

x(q,r) = q + Z}(r)V$o(q) , (64) 

where only the growing mode solution D{t) of Eq.(34) has been considered, and we 
introduced the potential $o(q), such that Vq$o = —(^o/Dq, where is the standard (i.e. 
Euclidean) Laplacian in Lagrangian coordinates; more in general, at this perturbative 
order covariant and partial derivatives with respect to the coincide. The potential $o 
is easily related to the initial peculiar gravitational potential defined in the Introduction, 

$0 = -(47rGa2pob^o)"Vo- 

More interesting is to derive from the above expression the corresponding Zel'dovich 
metric. It reads 

7f j^^(q, r) = 6,s {S\ + D{t)%,\{cD^ (6% + D{T)%,%{ci)^ . (65) 

One can of course diagonalize this expression by going to the principal axes of the 
deformation tensor. Calling Aq, the eigenvalues of the matrix $0)"/35 one finds 

7r^(q,r) = [l + /^(r)A.(q)f . (66) 
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Note that, contrary to what has been commonly done so far in the hterature, the 
metric tensor must be evaluated at second order in the displacement vector, in order 
to obtain back the correct Zel'dovich expressions for the dynamical variables (density, 
shear, etc ...). 

The above diagonal form of the metric allows a straightforward calculation of all 
the relevant quantities. The well-known expression for the mass density is consistently 
recovered, 

0^'''' = QbUi^ + DX^)-\ (67) 

a 

The peculiar velocity-gradient tensor has the same eigenframe of the metric; its 
eigenvalues read 

By summing over a the latter expression we can obtain the peculiar volume-expansion 
scalar 



and then the shear eigenvalues 

- y """" . (70) 



ZEL D'Xn 1 x-^ D'Xr 



The (conformally rescaled) electric tide = a'^E°'^ comes out just proportional to the 
shear. Its eigenvalues read 

cZEL /l^/^^2^ ^ZEL /V-rN 

= -47rGa gi,— . (71) 

These expressions for the shear and the tide completely agree with those obtained by 
Kofman and Pogosyan |^T| and Hui and Bertschinger |T^. The fact that metric, shear 
and tide have simultaneous eigenvectors shows that fluid elements in the Zel'dovich 
approximation actually evolve in a way much similar to a "silent universe" 0,0, i.e. 
with no influence from the environment, except for that implicit in the self-consistency 
of the initial conditions. 



4. Beyond the Newtonian approximation 

Having examined all the aspects of the formalism in the Newtonian limit, we are now 
ready to proceed to the next perturbative order in 1/c^. The PN terms -^w^^i^^ 
in Eq.(37) should be thought as small perturbations superposed on a Newtonian 
background 7^/3. The fact that the three-metric in the Newtonian limit is that of 
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Euclidean space in time-dependent curvilinear coordinates q'", implies that we can apply 
most of the standard tools of linear perturbation theory in a flat spatial background 
(actually, in an Einstein-de Sitter universe). In particular, we can classify the PN 
metric perturbations as scalar, vector and tensor modes, as usual. 
We then write 

,.,iPN) _ (PN)- , APN) l/.(PiV) APNh , (PN) .„21 

with 

^(PA^r^^^.CPA^r ^.(Piv^i^^O, (73) 

where greek indices after a single vertical bar, or nabla operators with a greek index, 
denote covariant differentiation in the Newtonian background three-space with metric 
7q,/3. In the above decomposition x^^^^ ^'^d (^(^^) represent PN scalar modes, ^^'^^ 
PN vector modes and vr^^^"* PN tensor ones (indices being raised by the contravariant 
background three-metric). We deliberately used the same symbols as in Section 2, in 
order to emphasize the analogy with the linear problem. Some of these PN modes, 
namely x*^^'^-' and 7r^^^\ also have a non- vanishing linear counterpart, as noticed in 
Section 2 (actually the linear part of tt^^^'' appears as a gauge mode in the equations), 
while others, namely and ^^^"^ are intrinsically non-linear. Unlike linear 

perturbation theory in a FRW background, metric perturbations of different rank do 
not decouple: this is because our time-dependent Newtonian background enters the 
equations not only through the metric 'jap, but also through the peculiar velocity- 
gradient tensor -(9"^, which also contains scalar, vector and tensor modes. This fact 
leads to non-linear scalar-vector, scalar-tensor and vector-tensor mode mixing, which 
also explains why we had to account for the vector modes ^j^''^'^ in the expansion of 
'w']^l^\ in spite of the irrotational character of our fluid motions. 

As in every perturbative calculation, some of the equations have the property to 
mix different perturbative orders. This is of course necessary in order to make the n- 
th order coefficients of the expansion calculable in terms of those of order n — 1. In 
our case the energy constraint and the extrinsic curvature evolution equation (which at 
the Newtonian level implies 7^°^ (7) = 0) play this role. Therefore we assume that the 
Newtonian metric and its derivatives arc known by solving the Raychaudhuri equation 
and the momentum constraint, and we calculate the PN metric perturbations in terms 
of them. 

Let us first compute the tensor V"^ = 47^"^ — (7^ + 2fi;)5"^ to first order in 1/c^. We 
obtain 

c^V%w^^^^) = -2(VV(^^)", + xP;) , (74) 
where now V^(-) = {■)\"'^- 
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In the energy constraint only the scalar x enters, 

V\ =\[^'- + - ^^Ga'e.iT''' - 1) , (75) 

where, here and from now on, we have dropped the superscript (PN) on PN terms. One 
can also obtain an equation for x from the trace of the evolution equation, Eq.(17), which 
is however equivalent to the latter, thanks to the Newtonian Raychaudhuri equation. 

The tensor perturbations vr"^ are instead determined via the evolution equation, 
Eq.(17) (actually from its trace-free part), 

V V", = 2^-; + 2 (2^ + ^) - \5^, - r,r^) -x\%. (76) 

A by-product of the latter equation is that linear tensor modes, which in the c ^ 00 
limit appear as harmonic functions (i.e. pure gauge modes), do not contribute to the 
r.h.s., i.e. to the Newtonian evolution of the system, as expected. 

In order to get an equation for the tensor modes decoupled from the scalar mode x 
we can resort to the equations obtained at the end of Section 3. To this aim we define 
the auxiliary Newtonian potential through 

V2^W = -1(^2_^M^^.^^ . (77) 

Using this definition in Eq.(75), we obtain 

X = -^f^ + - 2^(^) (78) 

a " 

(which, at the linear level reduces to the expression of Section 2, xo — —■^'^o)- Using 
this expression in Eq.(76) and replacing 1?"^ from Eq.(62) we obtain 

V^tt"^ = V" + 5"^ V^) *f ) + 2 (^^^"^ - ^''^^^^ , (79) 

which has the significant advantage of being explicitly second order (in any possible 
perturbative approach). This equation is one of the most important results of this 
paper: it gives (in the so-called near zone) the amount of gravitational waves emitted 
by non-linear cosmological perturbations, evolved within Newtonian gravity. In other 
terms, this equation, which is only applicable on scales well inside the horizon, describes 
gravitational waves produced by an inhomogeneous Newtonian background. 

We have '(9"^(q, r) = i^'(r)$o,"^(q), with -D(r) the growing mode solution of Eq.(34) 
(Z)(r) oc a(r) oc in the Einstein-de Sitter case) and $0 defined in Section 3. Therefore 



(80) 



where the symbol indicates the standard (Euclidean) form of the Laplacian in 
Lagrangian coordinates and ^0 = "^viTo) and indices are raised by the Kronecker symbol. 
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To completely determine the PN metric perturbations we still need the scalar mode 
( and the vector modes C,a, which can be computed through the momentum constraint. 
We then need -i?"^ at the PN order; it reads 

= + ^ (n^^, - + l^v 

By replacing this expression into the momentum constraint, and using the Newtonian 
identity Tup = ^^u\pi obtain 



-l^,a^,r + ^'/3i«n7 - n^/^r« + nv'a = , (82) 



which shows that ( and are implicitly determined by the Newtonian quantities, once 
X and vr"^ have been computed. 

5. Fluid— flow approach 

Following the fluid-flow approach, described in the classical review by Ellis (see also 
Ref . p3| ) , we can alternatively describe our system in terms of fluid properties, in our case 
matter density, volume-expansion scalar and shear tensor, and two geometric tensors, 
the electric and magnetic parts of the Weyl tensor deflned above. The derivation of the 
equations reported below is thoroughly described by Ellis [0 and will not be reported 
here. 

For most cosmological purposes it is convenient to adopt the conformal rescaling 
and FRW background subtraction described in Section 2. Therefore, we can start by 
writing the continuity equation directly in terms of the density contrast 6, 

^ + {1 + 6)^ = 0, (83) 

with denoting convective differentiation with respect to the conformal time r. In 
the Lagrangian frame, however, and for a scalar fleld, convective differentiation and 
partial differentiation coincide. The peculiar volume-expansion scalar obeys the 
Raychaudhuri equation which we can rewrite in the form 

r)i9 n' 1 

^ + + ^2 ^ B ^ 4^c'a2f?,5 = , (84) 
Dt a S 

where cr"^ = ^9°g — ^S'^p^ is the shear tensor. The shear, in turn, evolves according to 
^ + ^^"^ + l^^", + ^"y, - + £", = 0, (85) 
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where we have rescaled the electric tide as E'^^ = a'^E'^^, which can be written in terms 
of our new variables as 

= + ^^^"/^ + ^^"/^ - <yv\ + (n - ^^"/3^) • (86) 

Note that, for a generic second rank tensor A"^, one has 

- + 7^ /3 - ^ 7 ' (87j 

where ^ denotes the total derivative with respect to r, which in comoving coordinates 
coincides with the partial one. The two last terms in the r.h.s. come from writing the 
Christoffel symbols in our gauge. It is then clear that when the operator acts on 
either the shear or the complete t?"^ tensor, the second and third term in the r.h.s. cancel 
each other and the convective and total differentiation coincide. This cancellation also 
occurs for a generic A^^ if either the relevant Christoffel symbols vanish (as it is the 
case for the Newtonian limit in Eulerian coordinates) or the convective derivative acts 
on the eigenvalues of and such a tensor has the same eigenvectors of a"^ [as it is 
the case for the electric tide in the "silent universe" case [@|,|^. 
The electric tidal tensor evolves according to 

-\lP,{n''''n%\\s + r^'^\\\) + 47rGa2f.,(l + 5)a^p = , (88) 

where we have rescaled the magnetic tide as 7i"^ = a'^H'^^ and redefined the Levi-Civita 
tensor so that i]"'^"' = ^~i/2g"/37 (^for simplicity we used the same symbol after rescaling). 
Finally, the magnetic Weyl tensor evolves according to 

+1^^, (^""'^^ II. + ^""'^".ii.) = . (89) 

In the fluid-flow approach, besides the evolution equations, one has to satisfy a 
number of constraint equations. One has: the momentum constraint, which we rewrite 
in the form 

2 

-1 
3 

the Ti-a constraint (which we already used to define the magnetic tide in terms of 
derivatives of the spatial metric). 



^"/3||a = , (90) 
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the div £ constraint, 

= -lp,l.uV'^^<y\n", + ^-^a'Q,5,p , (92) 
and the div 7i constraint 

c'n^mc^ = i^,ia.v'''cr\£", ■ (93) 
The non-vanishing of div Ti, leading to the non-commutation of shear and tide, for 



the generic case of irrotational dust, follows from the recent analysis in Ref.[0, where 
it was shown that requiring div 7i = implies that Ti itself vanishes. 

In the above equations one also needs to know the three-metric ja/B- This can be 
obtained from the evolution equation 

7af3' = 27a7^^^ , (94) 

which is however only valid in Lagrangian coordinates. In order to completely fix the 
spatial dependence of the metric one also needs to specify the energy constraint (the 
trace of the Gauss-Codacci relations), which we rewrite in the form 

2 a' 

c\n - 6k) = a%af^^ - -7?^ - 4-^9 + IGTiGa^QbS . (95) 

3 a 

Although we will not use the fluid-flow approach in this paper it is interesting to 
have the complete form of the equations, with the correct powers of included, in order 
to understand the Newtonian meaning of the electric and magnetic tide. 



We are now ready to discuss the fluid-flow approach within the Newtonian 
approximation. We just have to discuss the order in our 1/c^ expansion at which 
the various tensors enter the equations above. It is immediately clear that the mass 
continuity equation, the Raychaudhuri equation and the shear evolution equation, where 
no explicit powers of c appear, just keep their form, once the various tensors are replaced 
by their Newtonian counterparts. So, we have 

5' + {1 + 6)^ = 0, (96) 

^' + -^ + + cx^^a^, + 4:7rGa^ghS = , (97) 

and 

+ ^a", + ha'^^ + a^^a^ - \6'^,a\a\ + B^, = Q. (98) 

On the other hand, by its very definition, the electric tide contains a contribution 
coming from the PN terms x ci^nd vr"^ because of the spatial curvature terms. It is 
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however immediate to realize that, once the expressions of Section 4 for the PN tensors 
X and vr"^ are used, one recovers the simpler form, 

^"^=(v°V,-^r^V^)¥.f), (99) 

which, in Eulerian coordinates reduces to the standard form 

Sab = <fg,AB - ^^Ai?V^.V9g . (100) 

On the other hand, if we replace in Eq.(91), the Newtonian peculiar velocity- 
gradient tensor, we obtain the well-known result (e.g. Ref. |]^) that the magnetic tensor 
identically vanishes in the Newtonian limit. This can be very easily shown by either 
applying the formalism of Section 3, i.e. writing -(^^g through covariant derivatives of the 
Lagrangian velocity potential, or by writing the same tensor in terms of the Jacobian 
matrix of Section 3. The physics underlying this result is the conformal flatness of the 
Newtonian spatial sections, implying the commutation of spatial covariant derivatives. 
A simple consequence of this fact is that, at the Newtonian level, the div S constraint 
reduces to 

^% = ^(^'QbS,^ , (101) 

which, owing to our expression for £, turns out to be just the gradient of the Lagrangian 
Poisson equation, namely 

VVf ^ = ^T^Ga^QbS . (102) 

Let us now come to the tide evolution equation. In that evolution equation the 
circulation of the magnetic tensor is multiplied by c^, which means that the PN part of 
curl Ti, is the source of non-locality in the Newtonian electric tide evolution equation. 
On the other hand, if we look at the magnetic tide evolution equation, which starts to 
be non-trivial at the PN order, we see that curl £ is consistently a PN quantity. 

The Newtonian meaning of the momentum constraint has been already discussed 
in Section 3. Also interesting is the div Ti. constraint, telling us that the general non- 
vanishing of the PN magnetic tensor (see also Ref.|[T^), implies that the Newtonian 
shear and electric tide do not commute, i.e. they have different eigenvectors (viceversa, 
their non-alinement causes a non-zero div 7-^). Another possible version of this result is 
that the ratio of the velocity potential to the gravitational potential beyond the linear 
regime becomes space-dependent. 

To summarize, we can say that within the Newtonian approximation the fluid- 
flow approach in Lagrangian coordinates can be formulated in terms of mass continuity, 
Raychaudhuri and shear evolution equations plus the Newtonian div £ constraint, which 
closes the system, provided we remind the circulation-free character of the electric tide 
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in this limit. Of course the direct use of a constraint to close the system of evolution 
equations, has the disadvantage of breaking the intrinsic hyperbolicity of the GR set of 
evolution equations, so that the entire method looses its basic feature. No ways out: 
this is the price to pay to the intrinsic non-causality of the Newtonian theory. 

The above discussion on the role of the PN magnetic tidal tensor, as causing 
non-locality in the Newtonian fluid-flow evolution equations, completely agrees with 
a similar analysis in Ref.[^. The main difference being that the results reported here 
have been obtained directly in Lagrangian space, while they worked in non-comoving 
(i.e. Eulerian) coordinates. A different point of view on the subject is expressed by 
Bertschinger and Hamilton |3^, according to which the magnetic part of the Weyl tensor 
is non- vanishing already at the Newtonian level. According to Ref . pl| , the difference 
might be "semantic" ; what is most important is that there is general agreement on the 
fundamental fact that the Newtonian tide evolution is affected by non-local terms. 

Finally, let me mention how the equations for the "silent universe" model 
are obtained in this frame. If in all the equations above we were allowed 
to set 7i"^ identically to zero, then all the spatial gradient terms would disappear from 
the evolution equations. In that case each fluid element evolves independently of the 
others, apart from the self-consistency of the initial data, and the whole dynamical 
problem is reduced to the solution of a set of nine (because the shear, the electric tide 
and the metric can be shown to be aligned [|1^) ordinary differential equations. The 
origin of this simplification is easy to understand from the physical point of view: setting 
the magnetic tide to zero forbids the presence of gravitational waves, which, given the 
absence of sound waves (because of the vanishing pressure), are the only way for our 
relativistic system to transport information among different fluid elements, during their 
evolution. However, as one could show by analysing the constraint equations above, the 
vanishing of the magnetic Weyl tensor requires very special initial conditions, which are 
incompatible with generic initial data for cosmological structure formation 0,0. On 
the other hand, on scales much larger than the horizon, the absence of causal correlation 
is precisely what we would expect to occur, so that the silent universe picture can be 
reasonably applied to describe the non-linear evolution of generic initial conditions on 
super-horizon scales (see Refs.||^,0]). In particular one can apply this model to study 
the asymptotic behaviour of large (super-horizon-size) patches of the universe in the 
presence of a non-zero cosmological constant; Bruni, Matarrese and Pantano showed 
that in such a model most of the universe volume expands in a de Sitter-like fashion, 
consistently with the Cosmic No-hair Theorem. Contracting fluid elements, which exist 
even in the presence of a cosmological constant, generically collapse to a triaxial spindle 
singularity (see also Ref. 0). 
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6. Conclusions 

In these notes we have reviewed a Lagrangian approach to the evolution of an irrotational 
and collisionless fluid in general relativity. The method is based on a standard 1/c 
expansion of the Einstein and continuity equations which leads to a purely Lagrangian 
derivation of the Newtonian approximation. One of the most important results in this 
respect is a simple expression for the Lagrangian metric in the Newtonian limit; this 
can be written in terms of the displacement vector S(q, r) = x(q, r) — q, from the 
Lagrangian coordinate q to the Eulerian one x of each fluid element, namely 



dq» A dq'' 



(103) 



The spatial metric is that of Euclidean space in time-dependent curvilinear coordinates, 
consistently with the intuitive notion of Lagrangian picture in the Newtonian limit. 

Next, we considered the post-Newtonian corrections to the metric and wrote 
equations for them. In particular, a simple and general equation for gravitational- 
wave emission from non-linear structures described through Newtonian gravity was 
given in Lagrangian coordinates, which is has the relevant feature of being non-local. 
A simple way to deal with the problem is to transform the equation in Eulerian form, 
where it is easier to deal with the Laplacian operator (which has there the standard 
Euclidean form), obtain the Eulerian gravitational- wave tensor ttab and then go back 
to the Lagrangian expression through 7rQ,^(q, r) = J'\j'^pTTAB{^i(l, t), t). One obtains 
the Eulerian expression 

VIttab = + ^ABVl¥f^ + 2 (^Mab - M^^b) , (104) 

with V^^'^^-' = — 1('*^^ — d'^^d^j^). This formula allow to calculate the amplitude of the 
gravitational-wave modes in terms of the velocity potential which in turn can be 
deduced from observational data on radial peculiar velocities of galaxies, applying the 
POTENT technique [|| . 



In the standard case, where the cosmological perturbations form a homogeneous and 
isotropic random field, we can obtain a heuristic perturbative estimate of their amplitude 
in terms of the rms density contrast and of the ratio of the typical perturbation scale A 
to the Hubble radius rn = cH~^. One simply has 

SLs(-)\ (105) 



'"""KthJ 



as it can be easily deduced from Eq.(80), specialized to an Einstein-de Sitter model. 
This effect gives rise to a stochastic background of gravitational waves which gets a 
non-negligible amplitude in the so-called extremely-low-frequency band, around 10^^^ 
- 10"^^ Hz. On much smaller scales, where the effect might be even more relevant. 
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pressure gradients and viscosity cannot be disregarded anymore and the entire formalism 
needs to be largely modified. One can roughly estimate that the present-day closure 
density of this gravitational-wave background is 

n,UX)-6tJ—y . (106) 

In standard scenarios for the formation of structure in the universe, the typical density 
contrast on scales 1-10 Mpc implies that Qg^ is about 10~^ - 10~^. One might speculate 
that such a background would give rise to secondary Cosmic Microwave Background 
anisotropies on intermediate angular scales through a sort of tensor Rees-Sciama effect 
(Matarrese and Mollerach, in preparation). Let us finally stress that the amplitude 
of this gravitational-wave contribution, ~ (5^(A/rj^)^, is an important counter-example 
to the rule given in the Introduction, according to which relativistic effects should be 
proportional to ^g/c^ ~ 5{\/rHY- This is, in my opinion, an encouraging result, which 
should show to the skeptical reader that the application of relativistic cosmology to 
small scales can be much more than an elegant exercise. 
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